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Abstract 

Clusters of galaxies are powerful cosmological probes, particularly if their masses can be de- 
termined. One possibility for mass determination is to study the cosmic microwave background 
(CMB) on small angular scales and observe deviations from a pure gradient due to lensing of mas- 
sive clusters. I show that, neglecting contamination, this technique has the power to determine 
cluster masses very accurately, in agreement with estimates by Seljak and Zaldarriaga (1999). 
However, the intrinsic small scale structure of the CMB significantly degrades this power. The 
resulting mass constraints are useless unless one imposes a prior on the concentration parameter 
c. With even a modest prior on c, an ambitious CMB experiment (0.5' resolution and l^iK per 
pixel) could determine masses of high redshift {z > 0.5) clusters with ~ 30% accuracy. 
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I. INTRODUCTION 



Clusters of galaxies are powerful cosmological probes. The abundance of clusters is very 
sensitive to the amplitude of fluctuations in the mass density field. For example, the present- 
day mass distribution of clusters offers perhaps the cleanest way to measure the normaliza- 
tion of the matter power spectrum . As we gain the ability 
to probe the cluster abundance at high redshifts, we can hope to measure the evolution of 
this normalization, an evolution which is sensitive to the dark energy and its equation of 



state 



m 



li. 



Because the mass function varies so rapidly with mass, accurate mass determination is 
critical if cluster constraints are to realize their potential. Finding clusters is not enough: 
even if a method finds all clusters and returns no false positives, it still does not necessarily 
have constraining power. It is quite possible that one technique will be optimal for finding 
clusters while another is a more powerful probe of cluster masses. For example, radio sur- 
veys which can measure the Sunyaev-Zel'dovich effect are very promising ways of detecting 
clusters even at high redshift [l5|. But they do not necessarily yield an accurate mass esti- 



mator. In contrast, weak gravitational lensing has drawbacks as a c^ 



3,0, 
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uster finder 
2J. 



Q,li3 



but 



is potentially a powerful way of measuring mass 

While lensing of background galaxies has been studied in depth over the last decade, 
lensing of the cosmic microwave background (CMB) has not been explored as carefully, 
mainly because the angular resolution needed to see the predicted signal is only now being 



attained or planned 
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26|. Here I discuss how accurately cluster masses can be 
determined with observations of small angle (~ 1') observations of the CMB. 

On these very small scales, to a first approximation the background primordial CMB 
is a gradient. The photons passing close to the center of the cluster get deflected so that 
they appear to be originating further away from the cluster center. On the cool side of the 
gradient, this means the temperature is slightly hotter than it would appear without lensing, 



and on the 
distinctive 



lot side slightly cooler. So, subtracting off the dipole, the lensing pattern is quite 



as shown in Fig. ^ 



Although there are many possible contaminants to this signal, its distinctive morpho. 



allows one to hope that the signal can be extracted from any sources of noise 
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Here I examine how accurately CMB lensing can determine the masses of clusters, accounting 
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FIG. 1: Idealized temperature pattern due to lensing of a cluster. Here, the background CMB 
is assumed to be a pure gradient with strength Ty = 13/xK arcmin"^ and this gradient has been 
removed. The cluster, at redshift 1 with a mass m2oo = 

lO^^/i~^M0, has an NFW profile. The 
area shown is within the virial radius of 3.3'. Jaggedness at the edge of the map is an indication 
of finite size (0.5') of the (noiseless) pixels. 

for the small scale structure present in the CMB, i.e., the fact that the background CMB is 
not a pure gradient. There are two sources of this small scale structure: one is the power 
which remains in the damping tail of the primordial {z ~ 1100) spectrum and the other is 
lensing by structures along the line of sight. Figure El gives an indication of the strength of 
the signal and the contributions from the CMB, which I will call CMB noise. The CMB 
noise dwarfs the signal on "large" scales. We will see that this significantly impairs the 
constraining power of CMB-cluster lensing. 

Section II describes the model for the mass density I will use throughout and the ensuing 
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FIG. 2: Anticipated signal from the lensing of a massive cluster and noise in Z-space. The cluster 
has mass 10^^/i~^Mq and is at redshift 1. The light CMB curve is without lensing along the line 
of sight; the heavy CMB line includes this effect. The instrumental noise corresponds to a 0.5' 
resolution experiment with 1/iK noise per pixel. 

constraints if CMB noise were absent. Section III shows how CMB noise pollutes the signal 
and which modes retain their power. Figures IHIITUI summarize the main results of the paper 
and can be understood without reading the text. Finally, in the conclusion, I emphasize 
some of the simplifications made throughout. 



II. CONSTRAINTS NEGLECTING CMB NOISE 



Numerical simulations 
profile. 
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find that clusters obey a Navarro-Frenk- White (NFW) 
A 



p{r) 



(1) 



r(cr + 

with two parameters the virial radius and concentration c. The virial radius is the radius 
within which the enclosed mass, m2oo; is 200 times the average mass in a critical density 
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universe in the same volume. Therefore, it is convenient to use m2oo instead of the virial 
radius as the second parameter, 

SOOtt , , , 

"^200 = -^Pcr{zK (2) 

where the evolving critical density Pct{z) is expressed in terms of the Hubble rate, 
3H'^{z)/8'n-G. The amplitude of the density is then 

4[ln(l + c)-c/(l + c)]- ^ ^ 

Throughout I will take c = 5 since this is a typical value emerging from simulations js^, 
although it does not reflect the changes of c with redshift or the (slight) changes with mass. 

Given this mass distribution, we can ask how the observed signal T in the CMB is affected 
by the lens. Lensing can be included by Taylor expanding the temperature field around the 
unlensed field T: 

m = T{e- 56) ^ T{e) -se-^ 

. ■ ^ + - se.„ ■ ^) (4) 

where 69 is the deflection angle due to lensing along the line of sight due to the cluster (h 



subscript) and unassociated large scale structure (iss subscript). Seljak and Zaldarriaga 
I28 I realized that on small scales the primoridial CMB has little structure so dT/dO in the 
first term can be set to a constant. Aligning the y-axis with this vector and calling the 
constant Ty leads to 

m ^ -Se^yMTy + (t{9) - ■ ^) . (5) 

To be clear, d9c\,y{0) is the ^/-component of the deflection angle due to the cluster at angular 
position 6, and Ty is the primordial gradient at the center of the cluster, the derivative of 
the temperature field with respect to 6y. The term in parentheses in Eq. is the lensed 
CMB temperature with mean zero and variance determined by the set of C^'s appropriate 
for the chosen cosmology. Note from Fig. El that lensing due to large scale structure gives 
additional power to the CMB at the scales of interest. 

The deflection angle due to the NFW profile is known to be js^ 

o9d[9) = -B-j-gidiOc r^} (6) 

CTy ds 
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where the distances are angular diameter distances (e.g. the angular diameter distance 
between the source and the lens dsi is {xs — Xl)/{^ + Zs) in a fiat universe where x is the 
comoving distance out to the relevant redshift). The function g is of order unity and peaks 
when its argument is equal to 1.3: 



1 J ln(x/2) + 



X < 1 



n(x) = - < ^^"^ f7) 

^ [ln(x/2) + "/^-™i'/") x>l' 

Figure ^ shows the temperature pattern on the sky - with the background gradient 
removed - due to lensing of an NFW cluster with mass IO^^/i-^Mq at redshift 1. Note the 
distinctive double lobe pattern on either side of the cluster. Here I have set the background 
gradient {Ty) to (Xgrad = 13/iK/arcmin, its rms value in a fiat cosmology with = 0.7, Qb = 
0.04, h = 0.7, cTg = 0.9, r = 0.1, and primordial spectral index n = 1. The hot and cold spots 
on either side of the cluster are then 9/iK an arcmin away from the cluster center, dropping 
off to 7/iK at the virial radius of 3.3'. 

The mild drop of the signal means that the outermost regions of the cluster contribute 
most to the mass constraints when only instrumental noise is considered. To see this and to 
contrast with the more realistic situation when other sources of noise are included later on, 
consider the Fisher matrix which determines the constraints on the parameters m2oo aiid c: 

F., - a^.. 1^^^^ (C j,^.^^- (8) 

Here a, (3 label the two parameters (Ai = lnm200; ^2 = ln(c)), while i,j, which are summed 
over, label the N pixels in the map. If the noise matrix is diagonal with uniform noise (j„ in 
all pixels, which might be true if only instrumental noise is considered, then the fractional 
constraint on the mass for example reduces to a single sum over all pixels: 

o^n i=i \ Cinm2oo / i=i 
That is, each pixel contributes an amount fu^i to the mass constraint. Figure El shows that 
although the largest contribution comes from the pixels a distance 0.4r„ from the cluster 
center, when we sum over all pixels in an annulus [oc / drrfii{r)], the annulus with the 
largest radius has the most constraining power. 

If this conclusion held up as more sources of noise were included, it would have dramatic 
implications for which clusters were best suited to be studied by CMB lensing. To see this, 
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FIG. 3: Contributions to the mass constraint as a function of distance from the cluster center. 
Although the largest individual contribution comes from pixels much closer to the center than the 
virial radius as indicated by the solid curve, the pixels at the edge contribute the most in total 
because there are more of them. This is in the unrealistic case where only uniform instrumental 
noise is considered. 

I will examine the constraints on the mass and concentration of a cluster by using all pixels 
out to the virial radius. Beyond the virial radius, the profile begins to deviate from an 
NFW profile. While we have no way of knowing in advance where the virial radius is, one 
can imagine fitting for all pixels within say 10' and getting non-negligible signal only from 
distances smaller than the virial radius. As we will see in §111, this discussion is academic 
because CMB noise renders pixels near the virial radius much less useful. 

With this strategy, we would constrain the mass of a cluster at lower redshift more tightly 
than one (of the same mass) at high redshift because the low redshift object subtends a larger 
fraction of the sky. We can see this in Figure 0] which shows the 1-a constraints on m2oo 
from an idealized CMB experiment with o"„ = 1/iK and pixel size A = 0.5'. The constraints 
in this unrealistic situation are astounding, percent level errors on the mass, even after 
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FIG. 4: Shaded contours show 1-a projections of SM/M of a cluster from a CMB experiment with 
noise per pixel of 1/iK and pixel size 0.5'. These constraints do not include CMB noise. Dashed 
lines indicate values for the cumulative number of clusters (total number above the given mass and 
redshift) expected in the standard ACDM cosmology (erg = 0.85). 

marginalizing over the concentration. Note the trend that for fixed mass, the constraints do 
indeed get stronger as the cluster moves to lower redshift. The dashed lines in the figure 
show the cumulative number of clusters expected at each mass and redshift. For example, 
the point z = 0.2,m2oo = 5 x lO^'^h'^MQ is intersected (roughly) by the contour labeled 
"10000." In this cosmology then, we expect 10^ clusters over the whole sky with mass greater 
than 5 x IO^^H'^Mq at redshifts greater than 0.2. (These contours were obtained using the 
Jenkins mass function Q].) 

III. CONSTRAINTS INCLUDING CMB NOISE 

To include noise from the CMB, I generate a set of C;'s for the chosen cosmology including 
lensing due to large scale structure using CMBFAST j^^l- The new covariance matrix is 
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then 



21 4- 1 

C = C^°^'' + ^ -^QPiicos %) (10) 



where 6ij is the angular distance between pixels i and j, and (7"°'^'' includes instrumental 
noise, assumed to be diagonal with variance a^. 

In the presence of this new source of noise, it is useful again to ask which modes contribute 
most to the mass constraint. Now though, since the covariance matrix is not diagonal, we 
need to generalize Eq. Q. If we diagonalize the full covariance matrix 

W = U'C-^U, (11) 

with the eigenvalues contained in the diagonal matrix W then the Fisher matrix of Eq. (jH)) 
can still be written as a sum of weights, but now the sum is not over pixels but over modes. 
That is, the contribution of a given mode to the mass constraint is 

2 



= Wia 



2 

grad 



dse, 



d In m2oo 



(12) 



The eigenvector associated with this mode is UjiT{9j). 

Fig. El illustrates how the signal in many modes is contaminated by CMB noise. Consider 
first the (red) circles which show /n for all the modes in the absence of CMB noise. The mode 
with the highest constraining power is depicted in the left panel of Fig. IHl it is essentially 
a gradient. Indeed, you can quickly estimate the mass constraint (with fixed concentration) 



from this mode: /n is of order 2000, leading to 5m2oo/'"^2oo ~ l/v2000 ~ a few percent, in 
agreement with the result for a cluster with m2oo = IO^^/i^^Mq at 2; = 1 shown in Fig. ^ 
The filled squares in Fig. El portray a much different picture when CMB noise is accounted 
for. Now there is only one mode with fu > 10; it is shown in the right panel of Fig. (HI It 
is sensitive to the difference in the cluster signal in the outer and inner parts of the cluster. 
This mode cannot be reproduced by the CMB because the CMB lacks the necessary small 
scale power, but the signal in this mode is also quite a bit smaller. Recall from Fig. ^ that 
the signal changes only mildly from its peak to its value at the virial radius. 

CMB-Cluster lensing then is sensitive not to the strength of the signal in both the hot and 
the cold spots, but rather to the structure of the signal in each. This changes the calculus of 
moving to lower redshift. All pixels on one side of the cluster no longer are weighted equally, 
so there is no longer an advantage in the larger clusters at low redshift. 
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FIG. 5: Constraining power of each of the top 10 modes when CMB noise is not included (red 
open circles) and when it is (black solid squares). The modes are linear combinations of the pixels. 
While of order ten modes are useful in the absence of CMB noise, only one mode (mode 6 here) 
contributes significantly to the mass constraints when CMB noise is taken into account. 

Figure [7| shows how the constraints change on an individual cluster when CMB noise 
is taken into account. The innermost ellipse shows the constraint on a massive cluster 
before accounting for CMB noise. While the concentration is not pinned down with very 
high accuracy, the mass is. When CMB noise is introduced, there is a degeneracy between 
concentration and mass. The inset in Fig. [7| gives a hint as to the origin of this degeneracy. 
The concentration parameter affects the deflection angle, and hence the signal, only close to 
the cluster center. Increasing the mass, on the other hand, increases the signal everywhere 
within the virial radius. Without CMB noise, all pixels at distances 6 > 0.56v therefore 
contribute to help pin down m2oo- The concentration is then determined by the innermost 
region, with its much smaller number of pixels; hence the concentration is known much less 
accurately. In the presence of CMB noise the situation changes: now the most important 
mode measures the difference between the signal in the inner and outer regions. Increasing 
the signal in the outer region by raising the mass can be offset by increasing the signal 
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FIG. 6: Left panel: The mode which is most constraining in the absence of CMB noise, but which 
gets annihilated when this noise is included. This is mode 10 in Fig. [51 Right panel: The one mode 
(mode 6 in Fig. which contributes significantly to the constraints when CMB noise is included. 

in the innermost region by raising the concentration. Therefore, only the sum of the two 
parameters is well constrained. 

The mass-concentration degeneracy can be broken if we use the knowledge obtained from 



simulations |3^ 



39| about the relation between the two. Imposing a prior on c with a 25% 
uncertainty leads to the intermediate ellipse in Fig. [7| The errors on the mass of this cluster 
then get reduced to useful levels. Even a more conservative prior on c of 50% leaves the mass 
errors quite tight: the factor of two increase in the c prior leads to only a ~ 40% increase in 
the mass error. 

Figure IHl shows the fractional uncertainty on cluster masses when all forms of noise consid- 
ered above are included. The errors are above 50% for clusters with mass below 10^^/i~^Mq 
even with the ambitious assumptions for noise and resolution. Note also that the uncertain- 
ties now get larger as the cluster moves to lower redshift. Although the cluster abundance 
depends heavily on as, so the numbers in Fig. |H1 (which assume as = 0.85) are only a guess, 
it is unlikely that there will be more than a handful of clusters for which CMB lensing will 
determine the mass to better than 50%. One possibility for a small scale CMB experiment is 
to sit on these handful (which are first identified in optical or Sunyaev-Zel'dovich surveys). 
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FIG. 7: Constraints on the mass and concentration of a 10^^/i~^Mq cluster at redshift 1 with pixel 
size 0.5' and noise l^K. Inner solid (red) ellipse does not account for CMB noise; largest (black) 
ellipse does; and medium-sized (blue) ellipse includes a 25% prior on concentration c. Inset in 
upper right demonstrates why CMB noise induces degeneracy between c and m2oo- Shown are the 
derivatives of the deflection angle with respect to these two parameters. Note that the deflection 
angle - and hence the signal - is sensitive to the concentration only close to the cluster center. 



On the other hand, CMB-cluster lensing does enable us to determine a combination of the 
mass and concentration. One way to exploit this is to impose a prior on the concentration. 
Fig. ini shows the resultant mass constraints with a 25% prior on the concentration. With 
this prior, we can obtain masses accurate to 50% for tens of thousands of clusters with 
masses above 2 x 10^^/i~^Mq at redshifts above 0.4. The uncertainties go down to less 
than 30% for the hundreds of clusters at redshifts greater than ~ 0.8 and masses above 
~ 5 X IO^^/i^^Mq. These constraints could prove very useful in determining dark energy 



parameters. Alternative 
context of a halo model 



one could use the constraint on the mass/concentration in the 
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FIG. 8: Fractional uncertainty on cluster mass from a CMB experiment with pixel size 0.5' and 
noise 1/iK. Contours are as described in Fig. |^ The constraints here include CMB noise. 

The above constraints were for an ambitious CMB experiment with noise per pixel of 
1/iK and 0.5' pixels. Figure ^1 shows how well less ambitious experiments would do in 
determining the mass of a lO^^/i^^M© cluster at z = 1. Apparently pixel sizes must be 
smaller than ~ 1.5' and noise per pixel smaller than 5/iK for any reasonable constraints 
on cluster masses even including a concentration prior. This would seem to rule out CMB 
lensing in Planck as a useful tool for measuring cluster masses. 



IV. CONCLUSION 



Since gravitational lensing is directly sensitive to the mass of a cluster, it seduces us 
into hoping that we can reject more traditional phenomenological mass indicators such as 
X-ray temperature and galaxy counts. In the case considered here, lensing of the CMB, 
unavoidable realities spoil some of this allure. Since the primordial CMB does have some 
small scale structure, this CMB noise significantly degrades the constraining potential of the 
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FIG. 9: Same as Figure |S1 except that here a 25% prior is imposed on concentration c. 

lensed field. The constraints may still prove useful in constraining cosmological parameters 
such as the dark energy equation of state, but we will have to include this degradation when 
evaluating different experimental proposals. 

While I believe this work is an important step on the road to assessing the utility of 
CMB-cluster lensing, it is incomplete in several ways. First, I have neglected a number of 
other contaminants such as point sources, the thermal Sunyaev-Zel'dovich effect, and the 
kinetic Zel'dovich effect. It is at least conceivable that these other effects differ sufficiently in 
frequency dependence and spatial morphology from the cluster lensing signal that they will 
not further degrade the projected constraints much. For example, if the morphology of the 
kinetic Sunyaev-Zel'dovich effect is much different from that depicted in the right panel of 
Fig. ini then it would have no impact on the mass determination, for we have seen that this 
is the mode which most tightly constraints the mass. Second, I have assumed a particularly 
simple parametrized form for the cluster mass distribution. Accounting for asphericity and 
substructure requires numerical simulations and these effects might add further 

noise to the mass estimators. Note that all of these effects would serve only to degrade the 
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FIG. 10: Constraints on the mass of a cluster at redshift 1 with m2oo = 10 /i Mpc including 
CMB noise and a prior on the concentration parameter. The point at bottom left here corresponds 
to the (1, 10^5) point in Fig. El 

mass constraints further. 

To close on a positive note, I point out that CMB experiments have consistently exceeded 
our expectations. There is thus every reason to believe that the high resolution, low noise 
experiments anticipated here will take place. If they do, then - along with input from 
numerical simulations such as the correlation between mass and concentration - we can 
hope to obtain interesting limits on the masses of thousands of clusters from CMB maps. 

Note: After this work was completed (but before it was circulated). Holder and 

n 
and Vale, Amblard, and White |30|] put out preprints studying very simi- 
lar issues. Although we reached the same general conclusion - realistic effects degrade the 
ability of CMB-cluster lensing to determine masses - our papers are by and large comple- 
mentary. While I have used analytic techniques and therefore simplified cluster models, the 
other two papers use more realistic numerical simulations. The simplicity has allowed me to 
provide quantitative measures of the efficiency of CMB-cluster lensing for a wide range of 
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cluster masses, redshifts, noise levels, and pixel sizes. Ultimately, these estimates will need 
to be reinforced with the realistic simulations performed in 0, [s^ . 
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